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Abstract 

We investigate the production process of the single scalar bottom quark 
(sbottom) and scalar charm quark (scharm) at HERA. The sbottom and 
scharm could be produced via an i?-parity breaking interactions A' 123 ^ 
in the Minimal Supersymmetric Standard Model (MSSM). These pro- 
cesses give a slight excess in the invariant mass M eq distribution of the 
high Q 2 deep inelastic scattering for sufficiently heavy gauginos and the 
scalar top (stop). For the light gauginos and stop, it is shown that the 
tagging of 6/c-quarks and charged leptons with high transverse energies 
will be indispensable to search for the processes. 
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1 Introduction 



The HERA accelerator is the only existing electron - proton collider with the center 
of mass energy of 300GeV (318GeV since 1998) produced by 27.5GeV electrons and 
820GeV (920GeV since 1998) protons |1[. The luminosity has steadily increased and 
the integrated luminosity of e + p collisions during 1994 - 1997 has reached 36.5pb _1 
and 47.7pb~ 1 for HI and ZEUS, respectively. One of the most attractive features of 
the HERA is to provide a unique possibility to search for new phenomena revealing 
the physics beyond the standard model (SM) through rare topologies in the final 
states characteristic to electron - proton collisions. 

In this paper we investigate a single scalar bottom quark (sbottom) and a scalar 
charm quark (scharm) production which is feasible at HERA ep — > bX(cX). Our 
theoretical framework is on the basis of the Minimal Supersymmetric (SUSY) Stan- 
dard Model (MSSM) with i?-parity breaking (RB) interactions 

In 1997 both the HI || and ZEUS || collaborations reported an event ex- 
cess in comparison with the SM expectations in the deep inelastic scattering (DIS) 
e + p —> e + X at large x and high Q 2 . The news very much excited high energy physics 
comminity. Various possibilities to understand the anomaly have extensively been 
examined by theoreticians since then 0, |^, ffl . We have also proposed an interpreta- 
tion of the anomalous event by the scalar top quark (stop) production in the SUSY 
models with RB interactions J7|, |8|]. Contrary to initial expectation, the novel fea- 
tures have gradually faded in the whole data sample with increasing experimental 
data. However, this fact shows that HERA could have the potentiality exploring 
physics beyond the SM. 

Here we are concerned with ep —> b/cX in contrast with ep — > t X in a previous 
work. The experimental signature of the present process has a large variety corre- 
sponding to decay modes of the sbottom/scharm. When decay modes b — > ec and 
c — > eb are dominant, the slight excess will be understood by the present scenario. 
On the other hand, if the gaugino modes dominate over previous ones, high Pt lep- 
ton(s) and 6/c-quarks would be typical signatures of our processes. In this respect, 
we are interested in the high Pt leptons observed by the HI ||. 

Throughout the present work whole calculations of decay widths and cross sec- 
tions have been performed by using the GRACE-SUSY system, an automatic compu- 
tation program for SUSY processes [|ItJ. While the GRACE-SUSY system is originally 
designed to treat such elementary sub-processes as e + e~, eq and qq collisions, we 
have recently succeeded in extending to ep and pp collisions. We use the extended 
new versions which includes an interface to the PDFLIB too. For the parton distri- 
bution function we have used CTEQ4M [|i~l"|] . 



*In this paper the term "electron" is used to describe generically electrons or positrons if not 
specified. 
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2 Models and constraints 



We start from the interaction Lagrangian based on the MSSM with an RB interac- 
tion 

L = \' ljk (uf L d k P L e - dkB^PiUj) + h.c, (1) 

where Pl,r read left and right handed chiral projection operators, respectively. The 
interaction Lagrangian (JJ) has been derived from the general RB superpotential 
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W H = XijkULjE^ + X^LiQjD-k + \% k U^D^D%, (2) 

where i, j and k are generation indices. The first two terms violate the lepton number 
L and the last term violates the baryon number B. Incorporating RB interactions 
into the MSSM we have a possibility to unveil yet unresolved problems as (i) the 
cosmic baryon number violation, (ii) the origin of the masses and the magnetic 
moments of neutrinos and (iii) some interesting rare processes induced by the L 
and/or B violation. The realization of the coupling among participating particles 
in Eq. ([!]) will be most suitable for the ep collider HERA because the squark ufi or 

dkR will be produced in the s-channel in e + -q sub-processes. 

e + + d k — > uf L (3) 
e+ + ii] -¥ d kR . (4) 

The upper bounds on the coupling constants A' ljfc have already been settled by, 
for instance, neutrinoless double beta decay II], charged current universality 



|i~2"| , |j~5l , atomic parity violation (APV) [[12], |T|, [T(| and v e mass |L7| . Some of possible 
nine coupling constants A' ljfc are severely constrained by experiments mentioned 
above. 

Here we pay attention to the single sbottom and scharm production, 

e + + c -> b R (5) 
e + + b -> c£, (6) 

which are realized via a non-zero -R-parity breaking coupling A' 123 . The most strin- 
gent upper bound on A' 123 comes from the charged current universality experiments 
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A' 123 ~0.1 (7) 



for m T ~ 200GeV. 



Sbottoms {hi, b R ) as well as stops are naturally mixed each other due to a large 
Yukawa coupling to their partner quark. The left- and right-handed sbottom {pL, 
b R ) are expressed in terms of the mass eigenstates (b\, b 2 ) through a mixing angle 

bh = \ cos 9b — b 2 sin 6 h (8) 
bn = bi sin 9b + b 2 cos 9b. (9) 
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We should note that the sbottom mixing is considered to be natural in the case 
of large tan (3, because the off-diagonal terms of the sbottom mass matrix are pro- 
portional to nib(Ab + /itan/3), where Af, and p, respectively, denote the trilinear 
coupling and the SUSY Higgs mass. Then the interaction Lagrangian (HD can easily 
be rewritten in terms of the mass eigenstates (61, 62)- in particular, it is worthy to 
note that both sbottoms b\ and b 2 could be produced in the eq collisions through 

h^ ec = — A' 123 (sin 6^61 e^P^c + cos Q^^P^c) + h.c. (10) 

We would emphasize that only one species of scharm c£ can couple to e + b in our 
scenario with A' 123 7^ 0. 

Interestingly, there exists a theoretical upper bound on the mass of the lighter 
sbottom bi. In the MSSM, since the left handed stop and sbottom bi form 
an SU(2) doublet, their masses include the same contribution from a soft scalar 
breaking mass mg 3 . We obtain an upper bound on the sbottom mass m| as 

m~ < m~ — m\ + ml — m^y cos 2(3. (11) 

The lighter sbottom b\ cannot be heavy for relatively light t 2 . As we consider the 
single sbottom and scharm production, we should take into account consistent model 
parameter sets allowing m^, m~ ~ 200GeV. For example, in order to obtain ~ 

200GeV, the heavier stop must be rather heavy ~ 250GeV. This means that 
the lighter stop must be lighter than about 250GeV. 

Next we examine the decay modes of the sbottom and scharm. We assume 
that the sbottom can decay through 61,2 — > ec, bx® and Wti, where Xi denotes 
neutralinos (i — 1 ~ 4). On the other hand, the scharm can decay into cl — > 
eb, cXi and sxt,2- The first decay mode for both sbottom and scharm occurs via 
the .R-parity breaking mode and the others proceed through i?-parity conserving 
interactions. 

At present, we know that the most stringent mass bound on the sbottom with the 
dominant decay mode 61 — > ec is > 242(204)GeV for Br(&i — > ec) = l (0.5). This 
has been obtained by the combined experimental results of the leptoquark searches 
by the CDF and DO experiments at Tevatron |18| . It should be noted that the mass 
bound could become lower when the branching ratio of the i?-parity breaking mode 
bi — > e c is smaller. 

We are aware of constraints from the precision measurements at LEP1. Poten- 
tially, contributions to the Ap from the ti and bi becomes large if their masses are 
not so different from the weak mass scale of rriz- We have checked contributions to 
the Ap from the ti and bj Jl9| could become small to such extent as 7 x 10 for 



mr = 200GeV, mr = lOOGeV and m T = 250GeV. 
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3 Numrical results 



In Table I we show two typical input parameter sets adopted throughout our numer- 
ical calculations. Output mass paramters of the neutralinos and the lighter chargino 
are also presented for reference. 

Table I Typical input parameter sets 



masses in GeV 


A 


B 


M 2 


300 


150 


tan (3 


12 


2 


fi 


-300 


-300 


m c L 


235 


230 


m t, 


250 


90 


m t~ 2 


400 


300 


0t 


0.42 


1.0 


Ob 


1.2 


1.0 


m~o 


143.4 


76.2 


Xi 






m~o 


254.4 


158.3 


X2 






m~ + 


254.4 


158.4 


XI 







For these parameter sets, we first calculate the decay branchig ratios of sbottoms 
and the left-handed scharm (see Table II), where we take A' 123 = 0.2[0.1] and = 
225[200]GeV for the set(A) [(B)]. 



Table II Branching ratios 





A B 


b\ — > e c 

&i -> wh 

bi total width 


0.688 0.084 
0.312 0.334 
0.0 0.181 
0.0 0.401 
0.090GeV 0.164GeV 


b 2 — > ec 
b 2 ^bx\ 

b 2 ->bxl 
b 2 -> Wti 
b 2 total width 


0.693 0.075 
0.307 0.095 
0.0 0.278 
0.0 0.551 
0.179GeV 0.352GeV 


c L — > e& 
cl -> cx? 
cl -> c x° 
cl -> s Xi 
cl total width 


0.959 0.006 
0.041 0.069 
0.0 0.291 
0.0 0.634 
0.191GeV 0.697 
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From Tables I and II, we find that the set (A) is characterized by the domi- 
nant decay modes &i j2 -> ec and cl e b. On the other hand, the SUSY decay 
modes bi^ — > bx® and b\ — > W t\ are dominated in the set(B). These properties 
are originated from the different values of the SU(2) gaugino SUSY breaking mass 
M 2 and the lighter stop mass m t ~. Large M 2 in the set (A) corresponds to heavy 
neutralinos \i ■ The i?-pariry conserving decay modes of the sbottoms/scharm are 
kinematically forbidden in this case. 



& 

o 

b - 0.1 



T 



T 



Q 2 > 10 4 GeV 2 
M eQ > 200GeV 



3L 
0.00 
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Figure 1: A' 123 dependence of the total cross section for e + p — > e + qX. We adopt 
parameter set (A) and = 225GeV. 

In Fig.l we show the A' 123 dependence of the total cross section for the process, 

e +p ^ e + g x, (12) 

where we take the set (A) and = 225GeV is assumed. To extract the relevant 
signal we adopt kinematical cuts Q 2 > 10 4 GeV 2 and M eq > 200GeV. In our calcula- 
tion two kinds of cuts on 9*, the angle between the outgoing and incoming positron 
in the eq rest frame, are introduced. As clearly seen from Fig.l, the cross section 
is more sensitive to the A' 123 for the more restrictive cut cos 9* < 0. In another 
words, more restrictive cut on cos#* would be efficient to extract the signal from 
background. 

In Fig. 2 we show the differential cross section against M eq for the set (A), A' 123 = 
0.2, m b ~ = 225GeV, Q 2 > 10 4 GeV 2 and cos^* < 0.5. It will be seen from Fig.2(a) 



6 



10 B = i — i i — r-| — i — i — i — r-p — i — i — r-p — i — i — r-p — i — i — i — pi — i — i — i — pi — n — 



10" 




JO" Li i i i I 



I .... I 



I . . i i I 



I ... .J 



200 210 220 230 240 250 260 270 



M eq (GeV) 




200 210 220 230 240 250 260 270 



M eq (GeV) 



Figure 2: M eq distribution for set(A), A' 123 = 0.2, m 6 ~ = 225GeV Q 2 > 10 4 GeV 2 and 
cos 9* < 0.5. (a) e + p — > e + cX (solid line) and e + p — > e + bX (dotted line) including 
contributions from 6 li2 and c^, respectively, (b) the signal (solid line) and the SM 
background (dotted line) for e + p — > e + qX. 
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that only the sbottom/scharm production contributes to the relevant signal. As 
the masses of two sbottoms and left-handed scharm are in the set (A) almost de- 
generate, three peaks corresponding to their masses appear in the range of M eq = 
220 ~ 240GeV. Summing up all sub-process contributions we obtain the observable 
differential cross section represented by Fig.2(b). We can clearly see a slight excess 
of the cross section in the range of M eq = 220 ~ 240GeV. Apparent peak structures, 
however, cannot be identified. 




L- - " " "" ' I I I I I 

100 120 140 160 180 200 

m sbottomi (GeV) 



Figure 3: The sbottom mass dependence of the sbottom total production cross 
sections with different final states. The set(B) and A' 123 = 0.1 are assumed. 

The mass dependence of the total production cross sections of the sbottom 
with different final states is shown in Fig.3, where the set(B) and A' 123 = 0.1 have 
been adopted. We should note that Br(bi — > ec) ~ 0.1 and even the light sbot- 
tom m~~100GeV is not excluded from the leptoquark searches at Tevatron. Then 
signatures relevant to the sbottom production have a large variety of characteristic 
features. Relatively light sbottom m^~150GeV dominantly decays into bxi and the 
heavier sbottom also decays into W t\ and bx®- As the lightest neutralino Xi can 
decay into the ordinary particles e ± cb or bsu e , the decay chains leading to relevant 
signatures are written as follows, 

b\X -> bxiX -»• b(e ± cb)X, b{bsv e )X (13) 

-> bx° 2 X -> b{tTxT)X -> b(£ + r{e ± cb))X, b{£ + r{bsv e ))X (14) 

-> bxlX -»• b(bbxi)X -»• b(bb(e ± cb))X, b(bb(bsu e ))X (15) 

-> WhX -> (£ut)(c^)X -> (^)(c(e ± c6))X, (^)(c(6si/ c ))X (16) 
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High Pt charged lepton(s) and/or 6/c-jet(s) are typical signatures of these processes. 
Especially, multi-charged leptons and 6/c-jets are almost the SM background free 
signals and they could effectively serve us in the sbottom search. 





Figure 4: Transverse momentum distributions of the final muon .Pt(/-0 ( a ) an d 
the c-quark Px(q) (b) for the process e~p — > bX — > W~tiX — > (fx~u ll )(cXi)X 
(solid line). The set(B), mr = 200GeV and A' 123 = 0.1 are as- 
sumed. Dotted and dashed lines, respectively, correspond to the backgrounds 
e~p -> e esc W~qX — > e esc {pTv^)qX and e~p -> e esc W^"cX ->• e esc ( / u"z/ M )cX. 

In Fig. 4 we show the differential cross sections for transverse momentum of the 
final muon Pt{h) and the c-quark Pt{q) for the process (16) as an example. Here 
we take the set(B), = 200GeV and A' 123 = 0.1. When we use the inclusive cross 
section of the muon to extract the signal, the single W production, ep — > e esc WqX — > 
Zesci^v^qX , is severe background, where e esc means the final electron escaped into 
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the beam pipe. As seen from Fig.4, the background cross section is much larger than 
the signal. However, if we can tag the c-quark, the background should be restricted 
to ep — > e esc WcX — >■ e esc (/iz/^)cX and its cross section is comparable to the signal. 
In fact we can expect the event excess in both Pr(fj) and -Pt(c) distributions. In 
addition to the muon and the c-quark, the information of the 6-quark and/or the 
high Pt electrons (positrons) coming from the neutralino Xi decays is available for 
us. All these informations make easier to extract the signals which we are searching 
for. 

4 Concluding remarks 

We have investigated a possible scenario to understand an excess of large x and 
high Q 2 events in the e + p — >■ e + X observed by the HI and ZEUS experiments at 
HERA. Our reasoning is based upon the resonance production of the sbottoms and 
the scharm with an i?-parity breaking interaction in the framework of the MSSM. 
We have focussed our attention upon its broad mass M eq distribution characteristic 
of the HERA events in addition to large x and high Q 2 . Assuming almost degenerate 
three mass eigenstates b\, b 2 and c£, we have simulated the mass distribution on the 
basis of our specific scenario. In the case of large M 2 and m 4 ~ the present scenario 
has its validity. If this is not the case, the sbottom decay dominantly into lighter 
sparticles. Then signals are characterized by high Pt charged lepton(s) and/or b/c- 
jet(s). Now planning vertex detectors at the HI and ZEUS provide us the efficient 
b/c tagging. Together with a luminosity upgrade bringing integrated luminosities, 
for instance, to the level of lfb -1 , it is expected that an event excess will be confirmed 
and our MSSM approach opens a new horizon. 
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